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MHD WAVE MODES RESOLVED IN EINE-SCALE 
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Within the last decade, due to significant improve¬ 
ments in the spatial and temporal resolution of chro¬ 
mospheric data, magnetohydrodynamic (MHD) wave 
studies in this fascinating region of the Sun’s atmo¬ 
sphere have risen to the forefront of solar physics 
research. In this Chapter we begin by reviewing 
the challenges and debates that have manifested in 
relation to MHD wave mode identification in fine- 
scale chromospheric magnetic structures, including 
spicules, fibrils and mottles. Next we go on to discuss 
how the process of accurately identifying MHD wave 


modes also has a crucial role to play in estimating 
their wave energy flux. This is of cardinal impor¬ 
tance for estimating what the possible contribution 
of MHD waves is to solar atmospheric heating. Fi¬ 
nally, we detail how such advances in chromospheric 
MHD wave studies have also allowed us, for the first 
time, to implement cutting-edge magnetoseismologi- 
cal techniques that provide new insight into the sub¬ 
resolution plasma structuring of the lower solar at¬ 
mosphere. 


1. INTRODUCTION 

Due to its complex and dynamic fine-scale struc¬ 
ture, the chromosphere is a particularly challeng¬ 
ing region of the Sun’s atmosphere to understand 
[see, e.g.. Judge, 2006]. It is now widely accepted 
that to model chromospheric dynamics, even on 
a magnetohydrodynamic (MHD) scale, while also 
calculating spectral line emission, one must realis¬ 
tically include the effects of partial ionization and 
radiative transfer in a multi-fluid plasma under 
non-LTE conditions [e.g., Hansteen et al., 2007]. 
Within the past decade there has been a concerted 
international effort to try and advance our under¬ 
standing of this tantalising layer of the solar atmo¬ 
sphere, which is thought to be a key part of solv¬ 
ing the solar atmospheric heating problem. There 
have certainly been major advances in chromo¬ 
spheric observations from high spatial and tem¬ 
poral resolution space-borne and ground-based in¬ 
struments. These include Hinode [Kosugi et ah, 
2007] launched in 2006 and the Rapid Oscilla¬ 
tions in the Solar Atmosphere [ROSA; Jess et ah, 
2010a] multi-wavelength camera system based at 
the Dunn Solar Telescope (DST), which became 
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operational in 2009. Furthermore, in 2013 the 
Interface Region Imaging Spectrograph [IRIS; De 
Pontieu et ah, 2014] was launched with the specific 
task of studying the previously little explored UV 
lines formed in the interface region between the 
chromosphere and corona in unprecedented detail. 

There now seems to be general agreement that 
the Sun’s magnetic field is primarily responsible for 
plasma heating in its atmosphere, from the pho¬ 
tosphere, up through the chromosphere, interface 
region and finally into the corona. However, a 
rigorous debate is still ongoing as to which par¬ 
ticular plasma processes are actually responsible. 
Historically, popular proposed mechanisms such 
as MHD wave dissipation or magnetic reconnec¬ 
tion have had little direct observational evidence 
to support them due to the small spatial scales 
involved and the limited resolution of past instru¬ 
mentation. For example, earlier space-borne EUV 
and X-ray wavelength instruments launched in the 
I990’s, with their limited spatial and/or temporal 
resolutions, e.g., the SOlar and Heliospheric Ob¬ 
servatory [SOHO; Domingo et ah, 1995] and the 
Transition Region and Coronal Explorer [TRACE; 
Handy et ah, 1999], could only detect less frequent 
high energy reconnection/wave events. However, 
it is now apparent that the majority of heating in 
the solar atmosphere must be taking place on small 
spatial scales, certainly less than ~ 100 km perpen¬ 
dicular to the magnetic field direction. As a result, 
the rarity of large-scale reconnection/wave events 
related to, e.g., flares and CMEs means that these 
do not play a key role. 
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Regarding the chromosphere, since it can be ob¬ 
served from ground-based telescopes with much 
higher spatial and temporal resolutions than are 
currently possible for the corona, it has presented 
an opportunity to study the small-scale and ever¬ 
present dynamics of thin (< 1000 km wide) chro¬ 
mospheric magnetic features such as spicules, fib¬ 
rils and mottles. These are particularly visi¬ 
ble in narrowband spectral hlters such as Ha 
and Ca ii H and K, which have aided tremen¬ 
dously in the identification of various MHD wave 
modes propagating along such hne-scale chromo¬ 
spheric magnetic structures. Prior to these mod¬ 
ern high-resolution chromospheric observations, an 
abundance of MHD waves were expected to ex¬ 
ist in the Sun’s lower atmosphere since it is in 
essence an elastic/compressible medium permeated 
by strong magnetic fields that are constantly being 
stressed and perturbed by the magneto-convective 
motions generated below (see Chapter 25 for an 
overview of photospheric wave modes). In agree¬ 
ment with these predictions it has now been con¬ 
firmed that there are indeed both Alfvenic and 
magneto-acoustic wave modes propagating along 
chromospheric waveguides at all times [e.g., De 
Pontieu et ah, 2007; He et ah, 2009a, b; Morton 
et ah, 2012a; Kuridze et ah, 2012, 2013, to name 
but a few]. 

A significant breakthrough now is that with con¬ 
temporary multi-instrumental studies we can see 
how small-scale disturbances generated in the pho¬ 
tosphere impact on the higher atmospheric layers. 
For example, in the particular case of ubiquitous 
small-scale (~ 1000 km diameter) photospheric 
vortical motions, Wedemeyer-Bohm et al. [2012] 
traced the resultant energy transfer of associated 
magnetic tornadoes up through the atmosphere 
by exploiting simultaneous photospheric, chromo¬ 
spheric and coronal observations using the ground- 
based CRisp Imaging Spectropolarimeter [CRISP; 
Scharmer et ah, 2008] at the Swedish Solar Tele¬ 
scope (SST) and also the Atmospheric Imaging As¬ 
sembly [AIA; Lemen et ah, 2012] on board the So¬ 
lar Dynamics Observatory (SDO). SDO, launched 
in 2010, is especially useful in detecting the coro¬ 
nal signatures of energy propagation from the lower 
atmosphere since it continually observes the full so¬ 
lar disc, sampling UV/EUV emission about every 
10 seconds. This now makes it much easier for 
observers to compare co-spatial/temporal coronal 
data to photospheric/chromospheric data gathered 
from specific, limited duration ground-based cam¬ 
paigns. 

With all these new possibilities, there has been a 
step-change in all-encompassing solar atmospheric 
MHD wave studies in the past decade. Now in 
simultaneous multi-wavelength observations it is 


possible to observe abundant chromospheric waves, 
along with their photospheric drivers and associ¬ 
ated coronal signatures [e.g., Morton et ah, 2014]. 
However, before we can start considering the en¬ 
ergetics of the waves and their possible contribu¬ 
tion to plasma heating, the very first step that 
must be taken is to accurately identify which wave 
modes are being observed. This is not a trivial 
task and has actually been the cause of much de¬ 
bate since the first wave interpretations of Hin- 
ode’s chromospheric data [see, e.g., Erdelyi &: Ee- 
dun, 2007; Van Doorsselaere et ah, 2008]. However, 
without this knowledge it not possible to quantify 
what the wave cut-off behaviour and most likely 
damping mechanisms could be. In fact, specific 
MHD wave modes can have different dominant (or 
competing) frequency-dependent damping mecha¬ 
nisms, e.g., thermal conduction, radiative cooling, 
kinematic viscosity, MHD radiation, resonant ab¬ 
sorption and phase mixing [see, e.g., Aschwanden, 
2004]. This can result in wildly varying damping 
rates for different MHD wave modes and, in turn, 
their ultimate effectiveness for atmospheric heat¬ 
ing. Hence, accurate quantification of MHD wave 
energetics must be founded on an precise identifi¬ 
cation of the actual wave mode (or combination of 
wave modes) being observed, as documented in the 
recent review by Jess et ah [2015]. 

Most recently, the benefits of accurate wave 
mode identification have been threefold; first they 
allow us to more precisely quantify the chromo¬ 
spheric energy flux associated with each mode 
[Morton et ah, 2012a; Van Doorsselaere et ah, 
2014]; second, the broadband frequency infor¬ 
mation of MHD waves allows us to study the 
frequency-dependent wave damping [Verth et ah, 
2010; Morton et ah, 2014]; thirdly, it allows a com¬ 
plimentary approach to understanding the fine- 
scale plasma structure of the chromosphere by im¬ 
plementing magnetoseismological techniques [Ee- 
dun et ah, 2011b; Verth et ah, 2011; Morton et ah, 
2012b; Kuridze et ah, 2013; Morton, 2014]. These 
are all crucial gains that were only made possi¬ 
ble in the last decade due to the launch of Hin- 
ode, combined with significant improvements in 
ground-based spectroscopic/imaging observations 
(ROSA/DST and CRISP/SST) and polarimetry, 
most notably with the Coronal Multi-channel Po- 
larimeter [CoMP; Tomczyk & McIntosh, 2009]. 
Now, with the launch of SDO and IRIS, there is 
an urgent need to combine the best cutting-edge 
observational and modelling studies so that for¬ 
ward leaps in solar physics can be facilitated. This 
will open a whole new era of studying the heat 
generated through hne-scale plasma dynamics in 
the solar atmosphere, not just through the study 
of waves, but also, e.g., instabilities and nanohare 
heating events. The following Sections 2 to 4 deal 
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Figure 1. Example of an MHD kink wave observed off-limb in a spicule using the Ca ll H filter of 
Hinode/SOT by De Pontieu et al. [2007]. The panels demonstrate the plane-of-sky projection of the 
spicule’s motion. The largest (left-hand) panel shows a time-distance diagram taken at a slice approxi¬ 
mately perpendicular to the spicule axis. The smaller panels reveal sequential snapshots of the spicule 
and indicate the perpendicular slice as a white line. Image reproduced from De Pontieu et al. [2007]. 


with MHD wave mode identification in chromo¬ 
spheric waveguides. Then, in Section 5 we go on 
to consider energy fiux estimates of particular wave 
modes. Finally, in Section 6 we also detail how con¬ 
temporary chromosopheric wave observations have 
been exploited for the purposes of advancing the 
field of solar atmospheric magnetoseismology. We 
would like to note that in solar physics literature 
the term “magnetoseismology” is also known as 
“MHD seismology”, or just simply “seismology”. 
However, for the purposes of consistency we will 
employ the term “magnetoseismology” throughout 
this Chapter. 

2. MHD KINK MODE IDENTIFICATION 

Although quasi-periodic line widths and Doppler 
velocities were detected in solar spicule data as 
far back as the 1960s [e.g., Nikol’Sky &: Sazanov, 
1967; Pasachoff et ah, 1968; Weart, 1970], the lim¬ 
ited spatial and temporal resolutions of this era 
prevented the identification of specific MHD wave 
modes. Spicules have actually been of much in¬ 
terest to solar physicists for much longer. In fact 
they were first reported in scientific literature as far 
back as Secchi [1877]. Nowadays they are known to 
be thin jets of plasma channelled by the magnetic 
field in the Sun’s lower atmosphere [see, e.g., the 
review by Zaqarashvili &: Erdelyi, 2009]. However, 
what causes their formation and what their con¬ 
tribution is to plasma heating (if any) is still the 


matter of fierce debate [see, e.g., De Pontieu et ah, 
2011; Klimchuk, 2012]. 

Spicules, which are predominantly rooted at net¬ 
work boundaries, are seen off-limb in a 2D pro¬ 
jection as a “thick forest” in chromospheric lines 
such as Ha and Ca ii H and K. The first claim 
of kink wave detection in spicules was by Kukhi- 
anidze et al. [2006] using Ha Doppler data from 
the coronagraph and universal spectrograph based 
at the Abastumani Astrophysical Observatory in 
Georgia. However, it took high resolution imaging 
data from the Solar Optical Telescope [SOT; Sue- 
matsu et ah, 2008; Tsuneta et ah, 2008b] onboard 
Hinode to observe this particular MHD wave mode 
more unambiguously [De Pontieu et ah, 2007]. Pe¬ 
riodic motions, perpendicular to the direction of 
the magnetic field, were detected in spicules us¬ 
ing the Ca II H filter of SOT (see Figure 1). This 
strongly suggested that the main restoring force 
for these waves was magnetic tension and led De 
Pontieu et ah [2007], and later [He et ah, 2009a], 
to simply interpret them as Alfven waves [Alfven, 
1942] with their phase speed, ca, governed by the 
well known relation, 


where B is the magnetic field strength, // is the 
magnetic permeability and p is the plasma den¬ 
sity. Theorists immediately started debating the 
validity of this interpretation [see, e.g., Erdelyi & 
Fedun, 2007; Van Doorsselaere et ah, 2008]. The 
base objection was that Alfven’s linear wave and 
planar geometry model assumed the plasma to be 
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completely homogeneous, and was therefore not ac¬ 
curate enough to predict the observed properties 
of waves travelling through the Sun’s inhomoge¬ 
neous and finely structured atmosphere. In fact, 
spicules have a finite width (diameter < 1000 km) 
and likely have a substantial variation in plasma 
density transverse to the direction of the magnetic 
field [e.g., Beckers, 1968], further fuelling the de¬ 
bate as to whether bulk Alfven waves were the 
correct interpretation. 

Recent 3D MHD radiative transfer simulations 
suggest that spicules could be formed by a lo¬ 
calised and enhanced Lorentz force at their base, 
which squeezes the chromospheric plasma in such 
a way that it is thrown up to lower coronal heights 
[Martmez-Sykora et ah, 2011, 2013]. Such cutting- 
edge numerical modelling supports the idea of 
spicules being overdense relative to the ambient 
plasma. If one assumes that on average spicules 
represent thin magnetically dominated filaments of 
plasma, with chromospheric densities and temper¬ 
atures that penetrate into the corona, then the ob¬ 
served transverse motions indeed have to be mod¬ 
elled as MHD waves propagating along an over- 
dense flux tube relative to the ambient plasma. 
Models of this type in flux tube, i.e., cylindrical, ge¬ 
ometry have been around since the late 1970s [e.g., 
Zaitsev &: Stepanov, 1975; Wentzel, 1979; Wilson, 
1979, 1980; Edwin &: Roberts, 1983]. Deriving 
the dispersion relations in such models requires 
the physical constraints that the total pressure 
perturbations and normal velocity components be 
continuous at the flux tube boundary. This was 
a worthwhile advance on Alfven’s simple model 
since it allowed for more realistic geometry and 
the possibility of both transverse magnetic field 
and plasma density inhomogeneities. The derived 
dispersion relations resulted in a much richer va¬ 
riety of MHD wave modes than was possible in 
Alfven’s more simple model. Assuming that the 
equilibrium plasma variation in the azimuthal di¬ 
rection of the flux tube is negligible, it allows for 
Fourier decomposition in that direction, and ulti¬ 
mately wave mode categorisation in terms of the 
integer azimuthal wave number, m. For a cylindri¬ 
cal flux tube, the lowest order azimuthal wavenum¬ 
ber {m = 0) results in two distinct decoupled MHD 
axisymmetric wave modes, i.e., the incompress¬ 
ible torsional Alfven [see, e.g., Hollweg, 1978] and 
the compressible sausage [see, e.g., Nakariakov et 
ah, 2003; Aschwanden et ah, 2004] modes. Their 
specific defining physical properties and their re¬ 
cent identification in chromospheric observations 
are discussed in Sections 3 & 4. 

After m = 0, the next integer azimuthal 

wavenumber is m = 1 and this is associated with 
the so-called kink mode, which is the particular 


MHD wave mode under discussion in this Section. 
As shown in Figure 2 (right), the key feature of 
this m = 1 mode is that it is the only value of 
m that produces a bulk transverse displacement of 
the flux tube. The m = 0 and all higher order 
(m > 2) fluting modes do not do this. Because of 
the plasma structuring the kink speed depends on 
both the internal magnetic field strength (Hj) and 
plasma density (pi), as well as the external mag¬ 
netic field strength (Hg) and plasma density (pe)- 
In the zero plasma-/? limit equilibrium demands 
Be = Bi. In this magnetically dominated plasma 
regime the kink speed, c^, is described in the thin 
tube (or long wavelength) limit as, 

where B = Be = Bi. Note that the value of 
the kink speed lies between that of the internal 
and external Alfven speeds. The kink mode is 
highly Alfvenic since its main restoring force is 
magnetic tension [see the discussion by Goossens 
et ah, 2009]. As discussed by Van Doorsselaere 
et al. [2008], it is also only weakly compressible 
in the long wavelength regime, hence it is un¬ 
likely that intensity perturbations due to compres¬ 
sion/rarefaction could be readily observed for such 
waves. Actually, De Pontieu et al. [2007] did not 
report any intensity perturbations concurrent with 
the transverse waves observed in spicules. There¬ 
fore, if present, they must be very small relative 
to the background intensity. So it is worthwhile to 
note that the absence of detectable intensity per¬ 
turbations cannot be used as an argument to dis¬ 
count the presence of kink waves in favour of bulk 
Alfven waves. Erdelyi & Fedun [2007] and Van 
Doorsselaere et al. [2008] also criticised the Alfven 
wave interpretation of De Pontieu et al. [2007] from 
the flux tube perspective. Since Alfven waves must 
be torsional in flux tube geometry, they would not 
display the bulk transverse motions observed in the 
imaging data of Figure 1, i.e., the flux tube would 
appear stationary. This was the main argument 
put forward to support the kink mode interpre¬ 
tation over the initial Alfven wave interpretation 
suggested by De Pontieu et al. [2007]. 

Regarding observed properties of kink waves in 
spicules, in a case study of 94 events, De Pon¬ 
tieu et al. [2007] reported periods of 100 — 500 s 
and transverse velocity amplitudes on the order 
of 10 — 25 kms“^. A subsequent, but more lim¬ 
ited case study by He et al. [2009a], also with Hin- 
ode/SOT data, reported the presence of higher fre¬ 
quency kink waves with periods as short as 40—50 s 
but with similar velocity amplitudes to that esti¬ 
mated by De Pontieu et al. [2007]. An advance 
on the work of De Pontieu et al. [2007] by He et 
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Figure 2. Cartoon illustrating the MHD sausage and kink waves in a magnetic flux tube. The sausage 
wave (left), which has m = 0, is characterized by an axi-symmetric contraction and expansion of the 
tube’s cross-section. This produces a periodic compression/rarefaction of both the plasma and magnetic 
field. The kink wave (right) has an azimuthal wavenumber m = 1. This is notable since it is the only 
value of m that causes a transverse displacement of the flux tube. In contrast to the sausage wave, the 
kink wave displacement/velocity field is not axi-symmetric about the flux tube axis. The red lines show 
the perturbed flux tube boundary and thick arrows show the corresponding displacement vectors. The 
thin arrows labelled B show the direction of the background magnetic field. Image from Morton et al. 
[2012a]. 


al. [2009a] was an attempt to actually measure the 
upward propagation speed of the waves, resulting 
in values between 60 — 150 kms“^ up to 7 Mm 
above the solar limb. A larger case study of 89 
spicules by Okamoto k, De Pontieu [2011] found 
a median period of 45 s, which is more consistent 
with the period estimates of He et al. [2009a] than 
De Pontieu et al. [2007]. Interestingly, Okamoto 
& De Pontieu [2011] also reported that 59% of the 
waves were propagating upwards, 21% propagat¬ 
ing downwards and 20% standing. Okamoto k De 
Pontieu [2011] also estimated propagation speeds 
of 164 — 267 kms“^ up to 7 Mm above the limb, 
but stated that inferred speeds greater than 1000 
kms“^ above 10 Mm could not be physical. A 
possible explanation of this is that spicule inten¬ 
sity in chromospheric lines is too diffuse at higher 
altitudes to be reliable enough for wave studies. 

Concerning another important class of chro¬ 
mospheric fine-scale magnetic structures, statis¬ 
tical studies of kink waves in fibrils have been 
made in a series of papers by Morton et al. 
[2012a, 2013, 2014], Fibrils are low-lying elon¬ 
gated structures, most clearly seen on disc, that 
span supergranular cells [Foukal, 1971; Zirin, 1972]. 
Like spicules, they too have narrow widths (< 1000 
km) and therefore need the highest resolution in¬ 


struments available to analyse their wave proper¬ 
ties. To this end, a range of statistical studies 
were performed by Morton et al. [2012a] to ex¬ 
ploit the fantastic capabilities of ROSA [Jess et 
ah, 2010a] equipped with a narrowband (0.25 A) 
Ha filter. Morton et al. [2012a] found that bulk 
transverse oscillations in fibrils, as in spicules, were 
omnipresent. It was reasoned that because fibrils 
appear dark in the line core of Ha, it was most 
likely a result of a density enhancement relative 
to the ambient plasma that causes the radiative 
emission from the photosphere to suffer increased 
dimming in their vicinity [Pietarila et ah, 2011]. 
Hence, Morton et ah [2012a] classed fibrils as over- 
dense waveguides and that their bulk transverse os¬ 
cillations must, along with spicules, be interpreted 
as kink waves, not Alfven waves. The combined 
studies of Morton et ah [2012a, 2013, 2014] anal¬ 
ysed 1688 fibril kink wave events in both the quiet 
Sun and active region chromospheres. The col¬ 
lated results gave most kink wave periods in the 
range of 94 — 130 s with velocity amplitudes of 
5 — 25 kms“^, which are certainly of the same order 
as that found in the aforementioned observations 
of kink waves in spicules. Perhaps this should be 
expected since both structures are rooted in inter- 
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TABLE 1. Average (or individually measured) properties of chromospheric MHD kink waves. 


Structure 

Region 

Max. displacement 

Period 

Max. velocity 

Kink speed 

No. Events 

Reference 



amplitude (km) 

(^) 

amplitude (kms”^) 

(kms-l) 



Spicule 

CH 

200 - 500 

150 - 350 

20 ± 5 


95 

De Pontieu et al. [2007] 


CH 

- 

60 - 240 

20 ± 5 

- 

- 

Suematsu et al. [2008] 


CH 

1000 

130 

15 

460 

1 

Kim et al. [2008] 



700 

180 

8 

310 

1 




800 

170 

9 

260 

1 



CH 

36 

48 

4.7 

75-150 

1 

He et al. [2009a] 



36 

37 

6.1 

59-117 

1 




130 

45 

18.1 

73 

1 




166 

50 

20.8 

109-145 

1 



CH 

55 ± 50 

45 ± 30 

7.4 ± 3.7 

160-305 

89 

Okamoto De Pontieu [2011] 



600 

180 

22 

- 

1 

Ebadi et al. [2012] 


QS 

670 

220 

19.2 

- 

1 

Jess et al. [2012] 



630 

139 

28.3 

- 

1 




160 

65 

14.8 

. 

1 




410 

158 

16.2 


1 




380 

129 

18.5 

. 

1 




200 

105 

11.8 


1 




190 

171 

7.2 

. 

1 



AR 

283 ± 218 

— 

14 ± 112 

- 

112 

Type-I -Pereira et al. [2012] 


AR 

463 ± 402 

— 

18 ± 12 

- 

58 

Type-II 


QS 

245 ± 211 

— 

16 ± 11 

- 

174 



CH 

342 ± 257 

- 

20 ± 12 

- 

170 


Fibrils 


135 

135 

1 

190 

1 

Pietarila et al. [2011] 


QS 

315 ± 130 

— 

6.4 ± 2.8 

50-90 

103 

Morton et al. [2012a] 


QS 

71 ± 37 

94 ± 61 

24.5 ± 1.8 

- 


Morton et al. [2013] 


QS 

94 ± 47 

116 ± 59 

5.5 ± 2.4 

- 

841 

Morton et al. [2014] 


AR 

73 ± 36 

130 ± 92 

4.4 ± 2.4 

- 

744 


Rapid Blue-shifted 


300 


8 


35 

Roupe van der Voort et al. [2009] 

Excursions 

CH 

200 

— 

4-5 

- 

960 

Sekse et al. [2012] 


QS 

200 

— 

8.5 

- 

1951 

average - Sekse et al. [2013] 



220 

— 

11.7 

- 

1951 

maximum 

Mottles 

QS 

200 ± 67 

165 ± 51 

8.0 ± 3.6 


42 

Kuridze et al. [2012] 


QS 

~ 172 

120 ± 10 

~ 9 

50 

1 

Kuridze et al. [2013] 


QS 

252 

180 ± 10 

8.8 ± 31 

101 ± 14 

1 



QS 

327 

180 ± 10 

11.4 ± 3.3 

79 ± 8 

1 



granular lanes where it is likely they are excited by 
similar/identical drivers. 

Kink waves on-disc have also been detected in 
mottles and via Rapid Blue-shifted Excursions ob¬ 
served in the blue wings of chromospheric spectral 
lines. Apart from their wave properties, mottles 
and Rapid Blue-shifted Excursions are of particu¬ 
lar interest to solar physicists since it is thought 
they may be related directly (or indirectly) to 
spicules. Mottels, like spicules and fibrils, are 
also thin magnetically aligned structures less than 
1000 km wide, but can either appear dark or bright 
in chromospheric lines. They were identified in 
scientific literature as far back as the early 1970s 
[e.g., Alissandrakis & Macris, 1971; Sawyer, 1972]. 
ROSA investigations of kink waves in mottles by 
Kuridze et al. [2012, 2013] revealed transverse ve¬ 
locity amplitudes on the order of 8 — 11 kms“^, 
which again are of the same order as found in 
spicules. The same is also true of kink waves de¬ 
tected in Rapid Blue-shifted Excursions in large 
scale statistical studies by Roupe van der Voort et 
al. [2009] and Sekse et al. [2012, 2013], where they 
found transverse velocity amplitudes in a similar 
range. The observational characteristics of kink 
waves in all the various chromospheric waveguides 
discussed in this Section are summarised in Ta¬ 
ble 1. 


3. MHD SAUSAGE MODE IDENTIFICATION 

Unlike the weakly compressible non-axisymmetric 
kink mode, the axi-symmetric sausage mode is 
highly compressible, producing periodic changes to 
the cross-sectional area of a magnetic flux tube, 
analogous to fluid motion driven in an elastic tube 
by a peristaltic pump, as shown in Figure 2 (left). 
The presence of such motion was first detected in 
the Sun’s lower atmosphere at the photospheric 
level in solar pores by Dorotovic et al. [2008] em¬ 
ploying the G-band filter of the SST, and subse¬ 
quently by Fujimura & Tsuneta [2009] using Hin- 
ode/SOT. These intense magnetic features are es¬ 
sentially like small sunspots without penumbrae. 
Employing ROSA G-band data, Morton et al. 
[2011] accurately measured the periodic area and 
intensity changes exhibited by solar pores. It was 
found that some pores exhibited a clear anti-phase 
behaviour between area and intensity oscillations 
that was strongly indicative of the sausage mode. 
For further and more in-depth discussions of photo- 
spheric sausage mode observations see Chapter 25. 
In the context of the current Chapter, these initial 
photospheric discoveries naturally led to the search 
for sausage modes higher up in the chromosphere. 
Using ROSA Ha data, Morton et al. [2012a] did in¬ 
deed detect the anti-phase behaviour of flux tube 
width changes and intensities in fibrils (see Fig¬ 
ure 3). Furthermore, these sausage waves were 
found to be concurrent with kink waves (whose 
identification in fibrils was previously discussed in 
Section 2). Morton et al. [2012a] found the sausage 
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Figure 3. Panel (a) depicts a cropped ROSA Ha snapshot containing a pair of dark, and hence dense, 
chromospheric flux tubes. Using the cross-cut (black line) to extract intensity information, panel (b) 
displays the resulting time-distance diagram revealing the dynamic motion of the waveguides. Times are 
given in seconds from the start of the data set, while the overplots are the results from a Gaussian fitting 
routine to show concurrent kink (red line shows the central axis of the structure) and sausage waves 
(yellow bars show the measured width of structure). Here there are counter-propagating kink waves with 
periods of 232 ± 8 s and phase speeds of 71 ± 22 kms“^ upwards and 87 ± 26 kms“^ downwards. The 
maximum transverse velocity amplitudes in both cases is about 5 kms“^. The sausage wave shown here 
has a period of 197 ± 8 s, a phase speed of 67 ± 15 kms“^ and maximum transverse velocity amplitude 
of 1 — 2 kms“^. Panel (c) displays a comparison between the detected intensity (blue) and width (red) 
perturbations resulting from the Gaussian fitting. Sausage waves can naturally cause such anti-phase 
behaviour. Image reproduced from Morton et al. [2012a]. 


waves to have periods in the range 135 — 241 s and 
transversal velocities on the order of 1 — 2 kms“^. 
The resulting comparisons that can be made be¬ 
tween kink/sausage energy fluxes in fibrils will be 
discussed in Section 5. 

Slightly earlier work by Jess et al. [2012], again 
with ROSA Ha data, actually detected similar 
joint kink and sausage mode signatures in spicules 
seen on-disc. However, Jess et al. [2012] did not 
explicitly associate the observed patterns as be¬ 
ing concurrent sausage/kink waves. The authors 
were primarily concerned with the footpoint driv¬ 
ing mechanisms of these chromospheric waves since 
their dataset had simultaneous photospheric (G- 
band) and lower-chromospheric (Ca ii K) image 
sequences. G-band photospheric intensity oscilla¬ 
tions showed a distinct phase difference of about 
90° across a magnetic bright point located at the 


footpoint of a group of chromospheric spicules. A 
2D MHD simulation was performed by Jess et 
al. [2012] using a compressive field-aligned foot¬ 
point driver of maximum amplitude 12.5 kms“^ 
combined with the observed spatial phase differ¬ 
ence. This actually resulted in combined kink and 
sausage waves similar to those observed in spicules 
by Jess et al. [2012] and fibrils by Morton et al. 
[2012a]. Hence, this could provide an explanation 
of why both highly and weakly compressive wave 
modes occur together in such hne-scale chromo¬ 
spheric structures. Thus far, sausage waves have 
not been detected in off-limb spicules. This could 
be due to complicated line-of-sight effects present 
in the “thick forest” of spicules seen off-limb. How¬ 
ever, for better understanding of spicule wave dy¬ 
namics and energetics, proving their existence (or 
not) should certainly be the focus of future studies. 
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Figure 4. Panels contain sequences of four CRISP Ha Doppler images of a spicule at different sam¬ 
pling times (in seconds, blueshift bright). The largest Doppler shifts in the spicule are approximately 
±55 km s“^. Image adapted from De Pontieu et al. [2012]. 


4. MHD TORSIONAL ALFVEN WAVE 
IDENTIFICATION 

Observations of small-scale photospheric vortical 
motions, detected via G-band bright point track¬ 
ing in intergranular lanes, have been the subject 
of much interest in recent years [e.g., Bonet et ah, 
2008, 2010; Wedemeyer-Bohm &: Rouppe van der 
Voort, 2009; Steiner et ah, 2010; Wedemeyer-Bohm 
et ah, 2012; Morton et ah, 2013]. G-band bright 
points, of the order 200 km in diameter, are of¬ 
ten co-spatial with kG magnetic flux concentra¬ 
tions. In such cases they are often referred to as 
magnetic bright points [e.g., Stenflo, 1985; Solanki, 
1993; Crockett et ah, 2009; Jess et ah, 2010b; Keys 
et ah, 2011]. Such small scale intense magnetic 
flux tubes rooted in vortex flow fields are natu¬ 
ral sources of torsional Alfven waves, as well as 
other MHD wave modes [e.g., Fedun et ah, 2011a; 
Shelyag et ah, 2013]. 

Torsional Alfven waves, if propagating in a near 
or sub-resolution flux tube could cause identifiable 
simultaneous periodic red and blue shifts in an ob¬ 
served spectral line. Since this process is not neces¬ 
sarily related to wave/energy dissipation inducing 
temperature fluctuations, it falls under the guise of 
periodic non-thermal spectral line broadening, pro¬ 
viding the torsional amplitudes are large enough 
to cause noticeable red/blue shifts [Zaqarashvili & 
Erdelyi, 2009]. With Ha data from the Solar Op¬ 
tical Universal Polarimeter (SOUP) based at the 
SST, Jess et ah [2009] detected such periodic spec¬ 
tral line broadening above a magnetic bright point 
group. Since there was an absence of both co- 
spatial intensity oscillations and bulk transverse 


motions, the periodic spectral line broadening was 
interpreted by Jess et ah [2009] as evidence of tor¬ 
sional Alfven waves. The estimated periods were in 
the range 126 — 700 s and the average line-of-sight 
velocity amplitude was approximately 1.5 kms“^. 

Theoretically, torsional Alfven waves can exist 
for any azimuthal wavenumber, m, and by defini¬ 
tion they are completely incompressible. A ques¬ 
tion that arises from interpreting observations as 
torsional Alfven waves is: How likely is it that 
such purely divergence-free MHD wave modes are 
actually excited in the photosphere? Follow up 
3D MHD numerical investigations by Fedun et 
ah [2011b] of flux tubes driven by vortex drivers 
demonstrated that although torsional Alfven waves 
could be the dominant wave mode, kink and 
sausage waves were still unavoidably present. This 
was due to the fact that the chosen spiral driver 
did not have the particular azimuthal symmetry 
(or asymmetry) of one distinct m value, but was 
in fact a superposition of different m values. Fea¬ 
ture tracking studies to estimate horizontal veloc¬ 
ity field components at vortex locations in the pho¬ 
tosphere [e.g., Morton et ah, 2013] show that this 
complicated scenario is actually much closer to re¬ 
ality. 

Since spicules are rooted in intergranular lanes 
where such vortex motion occurs, it is natural to 
assume that torsional motion should be common in 
all of these structures. To search for this, the high 
resolution capabilities of GRISP and the TRI-Port 
Polarimetric Echelle-Littrow (TRIPPEL) spectro¬ 
graph, both based at the SST, were exploited by 
De Pontieu et ah [2012]. They were successful in 
actually resolving red-blue Doppler velocity asym¬ 
metries across the width of spicules in both Ha 
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Figure 5. The left panel shows a cross-cut of the kink wave displacement field (arrows), perpendicular 
to the flux tube axis. The internal and external densities are discontinuous at the boundary, indicated 
by the circle. The right panel shows a kink “quasi-mode” displacement held (arrows) in a non-uniform 
tube with 1/R = 1. The non-uniform density layer between the internal and external plasma is indicated 
by the annulus. The resultant resonance will cause the transverse motion of the kink wave to be mode 
converted into localized non-axisymmteric (m = 1) torsional Alfven motion within the inhomogeneous 
layer, i.e., rotational motion will be amplihed as the displacement held evolves. Image adapted from 
Goossens et al. [2014]. 


and Ca ii H data, interpreting this as the clear 
signature of torsional Alfven waves (see Figure 4). 
Hence, a more complete picture is now emerging 
of spicule dynamics. These chromospheric and in¬ 
terface region magnetic structures support at least 
three distinct types of motion, i.e., field-aligned 
flows along with both kink and torsional waves 
[Sekse et ah, 2013]. Modelling all these motions 
as independent, De Pontieu et al. [2012] estimated 
that the best ht to observed data via Monte- 
Carlo simulations required field-aligned flows of 
50 — 100 kms“^, kink velocities of 15 — 20 kms“^ 
and torsional motions of 25 — 30 kms“^. A draw¬ 
back of this forward model is that by assuming 
these motions are independent, it completely ne¬ 
glects the physical magneto-fluid behaviour of the 
plasma. Realistically, these motions are coupled 
and this connectivity should be taken into account 
when interpreting data. 

In the particular case of the m = 1 kink wave, it 
was pointed out by Goossens et al. [2014] that the 
velocity field of this MHD wave mode, even with¬ 
out the presence of other modes, is actually a com¬ 
bination of both transverse and rotational motion. 
Although in the long wavelength limit the internal 
velocity field is purely transversal, as shown in Fig¬ 
ure 5 (left), the external field is dipolar in nature 
and could certainly contribute to rotational motion 
measured in observational data through artifacts 
of line-of-sight integration. Furthermore, this ro¬ 


tational motion can be significantly enhanced for 
the kink wave through the process of resonant ab¬ 
sorption shown in Figure 5 (right). In essence, 
this mechanism causes the transverse energy of a 
kink wave to be channelled to the m = 1 tor¬ 
sional Alfven wave in an inhomogeneous intermedi¬ 
ate layer between the internal and external plasma 
where the kink wave frequency matches the local 
Alfven frequency. 

Unlike the m = 0 torsional Alfven wave, the 
m = 1 torsional Alfven wave does not have az¬ 
imuthal symmetry [see, e.g.. Spruit, 1981] and 
would therefore not look the same to an observer 
from any given line-of-sight. Goossens et al. [2014] 
pointed out that the m = 1 rotational motions 
could produce very similar Doppler signatures to 
the m = 0 torsional Alfven wave if the observer’s 
line-of-sight is approximately perpendicular to the 
bulk transverse kink motion. Hence, if a spicule 
is observed to have a clear periodic transverse mo¬ 
tion, indicating the presence of a kink wave, the 
Doppler signal across its width is likely to have a 
significant contribution from its m = 1 rotational 
motion. Although De Pontieu et al. [2011] mainly 
interpreted red/blue Doppler asymmetries as be¬ 
ing due to the m = 0 torsional Alfven wave, they 
did not discount the possibility of m > 0 rota¬ 
tional motion being present. This offers a great 
opportunity for both theorists and observers to 
understand the interplay between different MHD 
wave modes and flows in fine-scale chromospheric 
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waveguides. Here, the transverse structure will be 
spatially resolved, with simultaneous imaging and 
Doppler data combined to give an accurate insight 
into the true nature of the plasma dynamics at 
work. 

So in summary, MHD wave mode identification 
in chromospheric waveguides has caused much and 
often heated debate. However, without this foun¬ 
dation knowledge we cannot truly understand the 
various important aspects related to their contribu¬ 
tion to the total energy budget of the solar atmo¬ 
sphere and their possible contribution to plasma 
heating. In this regard, the following Section 5 
reviews progress in determining more accurate en¬ 
ergy flux estimations of specific MHD wave modes 
observed in the chromosphere. 

5. MHD WAVE MODE ENERGY FLUX 

As discussed in Section 2, De Pontieu et 
al. [2007] interpreted the transverse motions of 
spicules not as MHD kink waves, but as bulk 
Alfven waves. This led them to also make energy 
flux estimates using the expression for such waves, 
i.e., 

E = ]^pv^CA , (3) 

where, v is the maximum transverse velocity am¬ 
plitude and the factor 1/2 comes from the time- 
averaged energy over one complete period. 

Equation (3) is only valid under the assump¬ 
tion of plasma homogeneity, which results in the 
equipartition between kinetic (KE) and magnetic 
(ME) energy. Eurthermore, De Pontieu et al. 
[2007] assumed that the energy associated with the 
transverse waves in spicules was of the same or¬ 
der elsewhere in the chromosphere and interface 
region, i.e., in regions where transverse waves were 
not observable due to the intensity S/N ratio be¬ 
ing too low. Hence, De Pontieu et al. [2007] as¬ 
sumed the filling factor of the wave flux energy to 
be unity, even though spicules themselves are esti¬ 
mated to have a filling factor of no more than about 
5% in the chromosphere [Makita, 2003; Klimchuk, 
2012]. This has led to a number of serious crit¬ 
icisms by various authors [e.g., Erdelyi & Fedun, 
2007; Van Doorsselaere et al., 2008, 2014; Goossens 
et ah, 2013]. 

For the reasons discussed in Section 2, it was 
pointed out by Erdelyi &: Fedun [2007] and Van 
Doorsselaere et al. [2008] that these transverse os¬ 
cillations are more accurately interpreted as MHD 
kink waves. Van Doorsselaere et al. [2008] also 
went on to point out that kink wave energy flux 
would be strongly influenced by the filling factor of 
the plasma structure it was propagating through. 
In fact, for an overdense flux tube relative to the 


ambient plasma, as can be assumed for spicules, 
the kink wave energy flux takes a maximum value 
within the tube itself and decays in the external 
region. The rate of decay depends on both the lon¬ 
gitudinal wavenumber and density contrast. The 
larger the longitudinal wavenumber or density con¬ 
trast, the faster the energy flux decays as a function 
of distance from the tube. 

Importantly, Goossens et al. [2013] demon¬ 
strated that for Pi ^ pe, a kink wave has no local 
equipartition of KE and ME. In the long wave¬ 
length approximation with Bi = Be, the ratio of 
ME to KE inside the tube is. 


and outside it is 


ME \ ^ Pi + pe 
KE J j 2pi 

ME \ _ Pi + pe 
KE J g 2pe 


(4) 

(5) 


Also from Goossens et al. [2013], the ratio of ex¬ 
ternal to internal total energy (TE = KE -|- ME) 
is, 

(TE)e _ 3pe + Pi /gs. 

(TE)j 2)Pi + Pe 

Equations (4) and (5) show that local energy 
equipartition is only possible if pe = Pi- An 
overdense flux tube (i.e., pi > pe) results in 
(ME/KE)j < 1 inside and (ME/KE)e > 1 out¬ 
side. Hence, KE dominates inside the tube and ME 
dominates outside. Note also that (TE)e/(TE)j < 
1, implying there is more of the total energy inside 
the tube than outside. 

To illustrate the spatial variance in the distribu¬ 
tion of energy with a specific numerical example, 
we take PijPe = 3. This results in (ME/KE)* ps 0.7 
and (ME/KE)e = 2. The ratio of total energies 
gives (TE)e/(TE)i = 0.6, hence, more than half 
the total energy is inside the tube. Denoting the 
flux tube radius as R, it can be shown that 90% of 
TE (= TEj -|- TEe) is within 2R of the flux tube 
axis, and 98.5% is within 577. Hence, even for a 
modest density ratio, pi/Pe = 3, this still results in 
a notable localised concentration of energy in the 
immediate neighbourhood of the flux tube. The 
larger the density ratio, the more localised this en¬ 
ergy concentration will be. Therefore, interpret¬ 
ing the transverse waves found in spicules as kink 
waves means we have to take account of the spa¬ 
tially varying nature of the energy flux. Van Doors¬ 
selaere et al. [2014] derived an expression for the 
spatially averaged kink wave energy flux in a multi¬ 
tube system based on the calculations presented by 
Goossens et al. [2013], assuming small filling fac¬ 
tors (/ < 0.1) as, 


E = + Pe)v^Ck 


(7) 
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where all the quantities are assumed to be aver¬ 
age values taken from kink waves propagating in a 
multi-flux tube system. In Equation (7) u is the av¬ 
erage kink wave maximum transverse velocity am¬ 
plitude at the various locations of the dense flux 
tubes. Since for a kink wave the transverse veloc¬ 
ity amplitude decays with distance from the flux 
tube, V in Equation (7) has a physically distinct be¬ 
haviour to the maximum transverse velocity ampli¬ 
tude shown in Equation (3) for bulk Alfven waves. 
In this simpler homogeneous plasma model the 
Alfven wave has a uniform velocity amplitude in 
space. Taking an upper bound spicule filling factor 
of / = 0.05, Van Doorsselaere et al. [2014] applied 
Equation (7) to the original bulk Alfven wave en¬ 
ergy flux estimates of De Pontieu et al. [2007] (de¬ 
rived using Equation (3)), and found they were re¬ 
duced from 4000—7000 Wm“^ to 200—700 Wm“^. 
This highlights the very important fact that if kink 
waves in an overdense solar waveguide are wrongly 
interpreted as bulk Alfven waves, it can lead to a 
substantial overestimation of the energy flux. In 
the particular case of De Pontieu et al. [2007], the 
overestimation is believed to be at least an order 
of magnitude. 

Regarding filling factors of fibrils seen on-disc, 
Morton et al. [2012a] estimated a comparable up¬ 
per bound to spicules of 4 — 5%. Por the energy 
flux estimate of kink waves in fibrils Morton et 
al. [2012a] only considered the energy inside the 
flux tubes, ignoring the external contribution. In 
essence this is similar to what was done by Van 
Doorsselaere et al. [2014] in deriving Equation (7). 
Interestingly, Morton et al. [2012a] estimated the 
kink wave energy flux to be 170 ± 110 Wm“^, the 
same order as that derived for spicules with the 
necessary filling factor correction by Van Doorsse¬ 
laere et al. [2014]. This should not be surprising 
since spicules and fibrils both have similar densi¬ 
ties, filling factors and transverse wave amplitudes. 

Incompressible wave energy flux in the different 
form of torsional Alfven waves above a magnetic 
bright point group was estimated by Jess et al. 
[2009] to be about 240 W m“^, assuming magnetic 
bright points cover at least 1.6% of the solar sur¬ 
face at any one time. This is of the same order 
as that estimated for filling factor corrected kink 
waves in fibrils and spicules. Again, this is not 
unexpected since Morton et al. [2013] showed that 
photospheric vortex motion, which is the natural 
driver of torsional Alfven waves, was also found to 
excite abundant chromospheric kink waves in Ha 
fibrils. 

Morton et al. [2012a] also estimated the sausage 
wave energy flux in fibrils to be on the order of 
460±150 W m“^, which is almost three times more 
than that found for kink waves, suggesting that 


compressive wave energy is more abundant than 
its incompressible counterpart. This has important 
implications for their ultimate fate, since compres¬ 
sive and incompressive MHD wave modes can have 
quite different physical damping mechanisms and 
rates, as discussed previously in Section 1. 

Here we add the caveat that such energy flux es¬ 
timates, filling factor arguments aside, are based 
on resolved wave amplitudes only. Therefore, a 
substantial amount of wave energy may still be 
unaccounted for. In the corona, it has already 
been suggested by many authors [e.g., Hassler et 
ah, 1990; Hassler & Moran, 1994; Banerjee et ah, 
1998, to name but a few], that measured non- 
thermal spectral line broadening between about 
20 — 50 kms“^ could have a significant contribution 
from sub-resolution waves. With respect to the 
chromosphere, the soon to be operational Daniel 
K. Inouye Solar Telescope (DKIST) in Maui, with 
its 4m diameter telescope, will offer a much im¬ 
proved tool to probe the smaller scale wave dynam¬ 
ics than is currently available (e.g., via the Im SST 
and the 0.76m DST). To put this into perspective, 
the 2-pixel diffraction-limited spatial resolution at 
Ca II K wavelengths obtained by the DKIST will 
be sub-20 km, compared with the near-100 km res¬ 
olution offered by the DST. By helping us to better 
understand the velocity fields and fine structure of 
chromosphere waves, it will also provide us with 
more accurate energy flux estimates. 

In summary, it was to be expected that from 
the preceding debate about chromospheric MHD 
wave mode identification, arguments would also 
arise about the actual energy flux they contain. In 
this Section we have highlighted some of the main 
differing ideas on this contentious issue. In the 
next Section, we go on to review how the discovery 
of these MHD wave modes has helped us advance 
the field of chromospheric magnetoseismology. 

6. ADVANCES IN CHROMOSPHERIC 
MAGNETOSEISMOLOGY 

Significantly, using Hinode/SOT Ca ii H data. 
He et al. [2009b] measured the variation in both 
propagation speed and velocity amplitude of kink 
waves as they travelled along spicules. Verth et al. 
[2011] exploited this detailed information for the 
purpose of implementing chromospheric magneto¬ 
seismology. Previously, magnetoseismology in the 
Sun’s atmosphere was limited to TRACE observa¬ 
tions of post-flare standing kink waves in coronal 
loops [e.g., Nakariakov & Ofman, 2001; Aschwan- 
den, 2004]. Since it was mostly the fundamental 
mode that was observed in such events, this did 
not provide enough information for wave theorists 
to determine how the plasma density and magnetic 
field were varying along such structures. Obtain- 
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Figure 6. The panels show magnetoseismically derived estimates of background variables along a spicule 
by Verth et al. [2011] in comparison with results of other authors. The left panel shows the estimated 
area expansion of the flux tube, normalized to unity at the spicule footpoint s = 0, where s is the held- 
aligned coordinate along the spicule. The dashed lines signify the 95% conhdence bounds. The dotted 
horizontal line denotes the upper limit estimated by Tsuneta et al. [2008a] using Hinode/SOT data. The 
middle panel indicates the estimated decrease in magnetic field strength with height, normalized to unity 
at s = 0. The dashed lines also signify the 95% confidence bounds. The dots show the average unsigned 
magnetic field strength from the radiative MHD simulations of De Pontieu et al. [2007] to a depth of 
1.4 Mm. The vertical dotted line indicates the position of the photosphere. The right panel illustrates 
both normalized plasma and electron density. The solid line indicates the magnetoseismically determined 
variation in plasma density from Verth et al. [2011], with the dashed lines indicating the 95% conhdence 
bounds. The crosses are from Table 1 of Makita [2003]. Regarding the normalized electron density, the 
circles are from Table XIX of Beckers [1968] and the triangles are from Table 1 of Makita [2003]. The 
estimate and uncertainties by Bjplseth [2008] are shown by the solid and dashed lines respectively. Image 
adapted from Verth et al. [2011]. 


ing more detailed information about field-aligned 
plasma inhomogeneity length scales from standing 
kink wave observations requires the detection of 
higher harmonics [see, e.g., Andries et ah, 2009], 
but unfortunately these were found much less fre¬ 
quently in the data. However, after the discovery 
of ubiquitous propagating kink waves in the chro¬ 
mosphere, this opened a whole new avenue in so¬ 
lar atmospheric magnetoseismology. The govern¬ 
ing kink wave equation that had initially been de¬ 
rived for standing kink waves in coronal loops of 
longitudinally varying magnetic field and plasma 
density [see, e.g., Ruderman et ah, 2008; Andries & 
Cally, 2011], could now be applied to observations 
of propagating waves in the lower atmosphere. The 
ordinary differential equation that describes the 
transverse velocity component of undamped kink 
waves in the thin tube regime is. 


d? / v\ 

ds^ Vi?/ c|(s) Vi?/ ’ 


( 8 ) 


where s is the magnetic field aligned co-ordinate, 
oj is the angular frequency, v{s) is the maximum 
transverse velocity component, i?(s) is the flux 
tube radius and, 




B\s) 


(9) 


where B{s) is the magnetic field strength, taken 
to be the same inside and outside the tube, and 
{p{s)) = [/9i(s) + /2e]/2 is the average of the inter¬ 
nal and external densities. 

If both the maximum transverse velocity, t(s), 
and kink speed, Cfc(s), are estimated from observa¬ 
tions, then Equation (6) can be solved for the only 
unknown, R{s). From the determined R{s) and 
the thin tube magnetic flux conservation relation, 

B(,) oc ^ . (10) 

the variation in magnetic field, H(s), along the flux 
tube can also be deduced. Combining B{s) with 
the original observational estimate of Ck{s), we can 
go back to Equation (9) for determining the field 
aligned variation in average plasma density, 

Ms)) oc ^ . (11) 

Therefore, observational estimates of v(s} and 
Cfc(s) allow us to determine the variation of both 
the magnetic field and plasma density along solar 
waveguides. 

Verth et al. [2011] pioneered this magnetoseis- 
mological approach to find the variation of R(s}, 
B(s) and {p{s)) along a spicule (see Figure 6). 
This technique was later implemented by Kuridze 
et al. [2013] and Morton [2014] in further investiga¬ 
tions of mottles and spicules, respectively. In fine- 
scale plasma structures of near-resolution width 
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such as spicules and mottles, R{s) can be diffi¬ 
cult to determine from intensity information alone 
[see, e.g., DeForest, 2007]. Also, traditional meth¬ 
ods for determining plasma density and magnetic 
field strengths in the chromosphere through spec¬ 
troscopy [e.g., Makita, 2003; Bjplseth, 2008] and 
polarimetry [e.g., Trujillo Bueno et ah, 2005; Cen¬ 
teno et ah, 2010] have their own inherent problems. 
Hence, magnetoseismology provides a much needed 
complementary approach in determining near- (or 
even sub-) resolution structuring of the chromo¬ 
sphere. 

Note that for Equation (6) the propagating wave 
envelope, v{s), is independent of uj since the effect 
of frequency-dependent damping is not included. 
Verth et al. [2011] pointed out that if damping was 
present this would result in underestimating the 
rate of change of R{s) and hence the other quan¬ 
tities, H(s) and (p(s)). In fact, the damping rate 
of kink waves and other MHD wave modes in the 
chromosphere are of much interest, but so far little 
is known. In contrast, the damping rates of post- 
flare/CME standing kink waves in the corona have 
been very well studied. In a statistical analysis 
of 52 standing kink wave events in coronal loops, 
using combined TRACE and SDO/AIA data, Ver- 
wichte et al. [2013] found most quality factors were 
in the range t/P ~ 1 — 4, where r is the damp¬ 
ing time and P the period. A widely supported 
physical mechanism to explain this is resonant ab¬ 
sorption [see the review by Goossens et ah, 2011]. 
In an overdense flux tube resonant behaviour is 
often modelled analytically by the inclusion of an 
annulus at the boundary where the value of den¬ 
sity decreases continuously from pi to pe- Hence, 
an Alfven continuum is introduced into the flux 
tube via the creation of a boundary layer. Since 
the kink frequency is between that of the inter¬ 
nal and external Alfven frequencies, at some posi¬ 
tion in the boundary layer the kink frequency will 
match that of the local Alfven frequency and a 
resonance will occur. This causes the kink wave to 
be mode converted to the m = 1 torsional Alfven 
wave in the boundary layer, resulting in the ob¬ 
served kink transverse motion becoming damped. 
Analytically, this process can be described most 
easily in the thin tube thin boundary (TTTB) ap¬ 
proximation, which predicts an exponential kink 
wave damping rate with a quality factor, 


L - p ^ Pi + Pe 

P I Pi- Pe' 


( 12 ) 


where I is the width of the inhomogeneous den¬ 
sity layer, R is the flux tube radius and the factor 
F depends on the functional form chosen for the 
decrease in density between pi and pe- To give a 
particular example, choosing a sinusoidal decrease 


results in E = 2/7r. It can be seen from Equation 
(12) that the damping time can be reduced by in¬ 
creasing the boundary layer width relative to the 
flux tube radius (larger l/R)-, and also by increas¬ 
ing the internal/external density contrast (larger 
Pi!pe)- The damping rate predicted by Equation 
(12) has been exploited to determine the cross-held 
variation in plasma density in solar atmospheric 
waveguides through observed damping rates. Pri¬ 
marily, this has been attempted for standing and 
propagating kink waves in the corona [see, e.g., As- 
chwanden et ah, 2003; Arregui et ah, 2007; Verth 
et ah, 2010]. Now there is such an extensive data 
set for the damping rates of coronal kink waves, 
even more advanced statistical models are now be¬ 
ing employed [e.g., Arregui et ah, 2013; Verwichte 
et ah, 2013; Arregui & Asensio Ramos, 2014]. In 
contrast, to date there have only been a few at¬ 
tempts at estimating the in situ damping rates of 
kink waves in the chromosphere [e.g., Kuridze et 
ah, 2012; Morton, 2014]. 

The basis of support for the mechanism of res¬ 
onant absorption to explain observed kink wave 
damping is mostly founded on two separate argu¬ 
ments. Eirstly, expected order-of-magnitude val¬ 
ues for both viscosity and resistivity in the corona 
would not account for the reasonably fast damp¬ 
ing rates [e.g., Nakariakov et ah, 1999]. Secondly, 
flux tubes in the solar atmosphere are unlikely to 
have perfect discontinuities in Alfven speed at their 
boundaries as idealised by Edwin & Roberts [1983]. 
Hence, inclusion of a more realistic Alfven con¬ 
tinuum between the internal and external plasma 
would naturally introduce a resonant layer for the 
kink wave. Eor a broadband frequency driver, reso¬ 
nant absorption and the resultant process of mode 
conversion will cause m = 1 torsional Alfven waves 
to be excited on many magnetic surfaces which 
will then phase mix. This may lead to Kelvin- 
Helmholtz instabilities between neighbouring mag¬ 
netic surfaces, which in turn will generate smaller 
length scales at which heating becomes more effi¬ 
cient [e.g., Ofman et ah, 1998; Terradas et ah, 2008; 
Antolin et ah, 2014]. 

Importantly, such broadband frequency propa¬ 
gating kink waves were discovered in coronal loops 
with the CoMP instrument by Tomczyk et ah 
[2007]. This inspired theorists to model the pro¬ 
cess of resonant absorption for propagating kink 
waves. Initial work in this area was by Terradas 
et ah [2010], who found that the damping length 
(Ld) for kink waves in the TTTB approximation 
is inversely proportional to the frequency, /, 

= y Ck . (13) 
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Hence, the process of resonant damping of kink 
waves acts like a low-pass filter in solar atmospheric 
waveguides. 

Verth et al. [2010] fitted Equation (13) to CoMP 
data of broadband frequency kink waves propagat¬ 
ing in coronal loops with ~ 600 kms“^ and 
found t/P ~ 2.7, consistent with the range ex¬ 
pected for standing kink waves {t/P ~ 1—4). Such 
frequency dependent damping should be detectable 
in velocity power spectra as a function of height in 
the solar atmosphere. In fact, Morton et al. [2014] 
searched for this by comparing velocity power spec¬ 
tra from both the chromosphere and corona using 
the ROSA/DST and CoMP instruments, respec¬ 
tively. Interpolating over the height spanned by 
the interface region between the ROSA/DST and 
CoMP data (approximately 10 — 15 Mm), it was 
found that the damping length of kink waves in the 
interface region was about 12.5% of that previously 
estimated in the corona. This provided tentative 
evidence of greatly enhanced damping of propa¬ 
gating kink waves in the Sun’s lower atmosphere. 
Morton et al. [2014] suggested that this could be 
caused by a combination of smaller quality factors 
and lower kink speeds in the interface region. The 
recently launched IRIS spacecraft, which is specifi¬ 
cally designed to give the highest spatial/temporal 
resolution yet in interface region spectral lines, will 
be an invaluable tool to push forward from these 
initial studies by Morton et al. [2014] and actually 
measure the changing properties of kink waves as 
they traverse this fascinating region. Certainly, a 
whole new era of chromospheric and interface re¬ 
gion magnetoseismology is opening up before us. 

7. SUMMARY 

Until the launch of Hinode in 2006, solar atmo¬ 
spheric MHD wave observers and theorists were al¬ 
most exclusively focused on the corona. However, 
with the discovery of ubiquitous transverse waves 
in chromospheric and interface region spicules by 
De Pontieu et al. [2007], this branch of solar at¬ 
mospheric wave research gained a whole new lease 
of life. Now there is vigorous debate about which 
particular MHD wave mode (or modes) are being 
observed in the Sun’s lower atmosphere. More¬ 
over, this has naturally lead to intense discussions 
about how best to quantify their associated en¬ 
ergy flux, knowledge of which is crucial for under¬ 
standing the contribution of waves to plasma heat¬ 
ing. The plethora of high spatial/temporal ground- 
and space-based chromospheric imaging and spec¬ 
troscopic data now available has also allowed the 
first magnetoseismic studies of this fascinating and 
complex region. This has proved an especially 
useful complementary tool to probe the fine-scale 


plasma structuring of the chromosphere. With the 
more traditional methods of spectroscopy and po- 
larimetry, it is very difficult to estimate even aver¬ 
age plasma densities and magnetic field strengths 
in near-resolution width and short-lived chromo¬ 
spheric features such as spicules, fibrils and mot¬ 
tles, far less how they vary in time and space. Now 
magnetoseismology really has something to offer 
in this regard. In the context of all-encompassing 
studies of MHD wave propagation and energy de¬ 
position throughout the whole solar atmosphere, it 
is becoming increasingly clear that chromospheric 
waves play a vital role. In conclusion, a whole new 
era of chromospheric MHD wave research has truly 
unfolded, offering fantastic opportunities and chal¬ 
lenges to both theorists and observers alike. 
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